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Abstract. General cosmological models with spinor and scalar fields playing the role of gravita- 
tional sources are analyzed. The Noether symmetry approach is taken as a criterion to constrain the 
undefined potentials and couplings of the generic actions. For all the found Noether symmetries the 
corresponding dynamical systems can be analytically integrated. The obtained cosmological solu- 
tions describe the early and late Universe as expected by basing on the known eras of the Universe. 
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INTRODUCTION 

In the inflationary paradigm [Ql 0, @] a scalar field called inflaton is the responsible for 
a rapid accelerated expansion of the primordial Universe. The inflation solves the prob- 
lems of flatness, isotropy of the microwave radiation background and unwanted relics. 
But the recent observations show that the present Universe also expands acceleratedly 
[HHH], which is not expected by the standard cosmology. The most popular attempt of 
explaining this acceleration is to suppose that a scalar field is also inflating the Universe 
today This strange fluid that inflates the late Universe is called dark energy. Ad- 

ditionally to this issue, there is an old problem related to the galaxy rotation curves 
the common matter cannot account for the observed dynamics of the galaxies and so 
another type of matter is needed D9|]. Its nature is yet unknown and several candidates 



are proposed in the literature [10]. The referred unknown matter is called dark matter. 

The possibility of a Dirac spinor field to be the source of the accelerated periods 
of the Universe has been explored in the literature! 1 1 , 12, OL [HL [01]. This kind of 



field also produces a standard matter behavior [|16l 11711 . which may work as a model of 
fermionic dark matter [fl8l fl9|] . In this scope we will analyze the general role played 
by spinor and scalar fields in Cosmology. For this task, it will be taken general actions 
without specifying the potentials and couplings of the fields a priori, since it would 
be generally done in an ad hoc way. Instead of this, the undefined functions will be 
constrained b y th e requirement that the action present a Noether symmetry. Such an 
approach [2(i 21, 22, 23, 24, 25] may work as an insight on the choice of the potentials 



and couplings for the models, at the same time that it is a useful tool for the integration 
of the field equations. 



The models will be analyzed for a flat Friedmann-Robertson-Walker metric with the 
signature (+, — , — , — ) and the associated scalar curvature defined by 



, a d 2 
R = 6[- + ^ 

a a L 



following the convention of the natural units %%G = c = fi= 1. 



AN OVERVIEW OF THE NOETHER SYMMETRY AND DIRAC 

SPINOR FIELDS 

The Noether symmetry approach 

Consider a Lagrangian of the form Jz? = Jzf(q k ,qk), which is associated with the 
energy function 

E* = (1) 

and the following point transformation 

cik — >Qk = qt + ea k (q), (2) 

where the q k = qk(t) are the generalized coordinates, the dot represents a time derivative, 
the a k are known functions of n variables - n is the number of generalized coordinates 
- and e is an arbitrary infinitesimal parameter. 

Any invertible transformation of the generalized coordinates like © induces a trans- 
formation of generalized velocities, such that 

G*(q) = ^-qi, (3) 
dqi 

where the matrix / = || d Q / dqi\\ is the Jacobian of the coordinate transformation, whose 
determinant is supposed to be non-null. The induced transformation © is represented 
by the vector field 

X ^( ff ^ft + dt dqt) ' ^ 

which is defined in the space tangent to the point of the transformation and is called 
infinitesimal generator of symmetry . The a k = CC fe (q) are the coefficients of the generator 
of symmetry. 

The Lagrangian Jzf is invariant under the transformation @ if 

L X ^ = 0, (5) 
with Lx designating the Lie derivative with respect to the vector field X. 



By using the E-L equations and © one can show the identity 



which implies that 

M o = £«a-3-, (7) 

where Mq stands for a conserved quantity (or constant of motion) associated with the 
transformation X. Then it follows that if the condition © is satisfied, the Lagrangian Jzf 



presents a Noether symmetry [|26l |27J 



It can be shown that the new set of variables {Qk(q)} obeys the following system of 
differential equations 



where = 1, 2, ... fc— 1 and the and z are the new coordinates Qk(t\), with z being a 
cyclic variable. 

The Noether symmetry approach consists in taking a given Lagrangian of the form 
Jzf = with undefined functions (generic potentials and/or couplings) and 

selecting their possible forms by the requirement that Jzf satisfy the condition ©, which 
is performed if there is a set of 0^ such that at least one of them is different from zero. 
So, if the functions that satisfy such a condition exist and are taken for the Lagrangian, 
it presents a Noether symmetry and there exists a conserved quantity - given by £7J) - 
associated with its dynamics and a cyclic variable is allowed - through equations ((8) -, 
which can help integrate the dynamical system. 



Dirac spinor fields in curved space-times 



The Lagrangian that describes a spinor field in a flat space-time 
time -, called Dirac Lagrangian, is the following 



Minkowski space- 



j [ ¥fda¥-{d a \j/)f\j/]-m\j/\j/-V, 



(9) 



where i// and xfr = i/z^y denote the spinor field and its adjoint, respectively - the dagger 
stands for the complex conjugate transpose -, m is the mass parameter and V is the self- 
interaction potential. When a coordinate transformation is performed, the referred spinor 
field transforms according to the Lorentz group |28|] in the following way 



y/(jj) = exp 



ab 



\j/(x) 



(10) 



where the X ab represent the Lorentz transformation parameters and the E are the 
generators of the spinor representation, which are defined by 

^ ab = \[f^\- (ID 
In this definition, the y" are matrixes 4x4 that satisfy the Clifford algebra 

{fy} = fyb + ybf =2T1 ab . (12) 

The generalization of the Lagrangian © for a spinor field in a curved space-time 
0113] reads 

^= 1 -[W^W-(D fl W)^V}-rnWY-V, (13) 
being the covariant derivatives defined by 

D^\j/=(d^-a^\ t f, (14) 

D^\j7 = d^W+¥^, (15) 

where £2^ is the spin connection 

= K x -e v b {d a e\) ] T°T\ (16) 

The generalized y matrixes are 

= ^f, (17) 

which satisfy the generalized Clifford algebra 

{r",r v } = 2^ v , (18) 

with being called tetrad or vierbein and defined in the form 

dB a 

g^e v b = r] ab , where e£ = — . (19) 

The tetrad is the object that carries the information of the generalized coordinate trans- 
formation to which the spinor field is subjected. 

According to the Pauli-Fierz theorem plQ . the potential V contained in (fT3l) is an 
exclusive function of the scalar bilinear *F = T/7yA and the pseudo-scalar T/7y 5 i//. By 
simplicity, in our models we will suppose that V is a function only of X V. 

The E-L equations for T/7 and y applied to (fT3l) generate 

dV 

lY^D^y-my- — = 0, (20) 
D At v7r^ t +mv7-^ = 0, (21) 

and so we have the Dirac equations for the spinor field and its adjoint in a curved space- 
time, i.e., in the presence of a gravitational field. 



SPINOR FIELD NON-MINIMALLY COUPLED WITH GRAVITY 



Point-like Lagrangian and field equations 

The general action of a model for a Dirac spinor field non-minimally coupled with the 
gravitational field through an arbitrary function of *F reads 

S = J V=gd 4 x{F(V)R + ^ [xjfT^D^-iD^xjf^Vf] ~V(W)}, (22) 

where R is the Ricci scalar. It will be clear later why it is not supposed a massive term 
in this action a priori. 

For a spatially flat F-R-W metric, one can obtain from (1221 - by means of an integra- 
tion by parts - the following point-like Lagrangian 

JSf = 6ad 2 F + 6a 2 cWF' + ^a 3 (ij/y yr - yrf yr) + a 3 V, (23) 

with the prime denoting a derivative with respect to x ¥. 

From the E-L equations for yr and yr applied to d23l . it follows the Dirac equations 
for the spinor field and its adjoint coupled with the gravitational field, namely, 

3 

yr + -Hyr + ifyrV' -i6(H + 2H 2 )fyrF' = 0, (24) 
3 

ij7 + -Hyr - ixjrfv' + i6(H + 2H 2 )yffF' = 0, (25) 

where H = d/a denotes the Hubble parameter. 

The acceleration equation follows from the E-L equations for a, yielding 



a _ pf + 3pf 



(26) 



a 12F 

where the expressions for the energy density pf and pressure pf are defined by 

p f = V- 6HF ,X ¥, (27) 
p f = [V' -6(H + 2H 2 )F , ] X V - V + 2(F'*i> + 2HF /X ¥ + F"^ 2 ). (28) 

By imposing that the energy function - see expression (0Q) - associated with (1231) 
vanishes, we obtain the modified Friedmann equation, i.e., 

Potentials and couplings from the Noether symmetry 



In terms of the spinor components, the Lagrangian (1231 can be written as 

4 4 

Jgf = 6ad 2 F + 6a 2 dF' £ e t ( yr? yr-, + yrf yr^j + l -a> £ ( yr* i# - yr* yr^j + a 3 V, (30) 



and the dynamics is now described by nine coordinates, being the configuration space 
of the system represented by (a, i//J, y/j, a, y/| , y/}), with j = 1, 2, 3, 4. 

Starting from the Noether symmetry condition, Lx«if = XJzf = 0, applied to d30l) . 
with X defined for the present dynamics as 




we arrive at an equation that depends explicitly on a 2 , dy*, dxffj, yjy£, ¥j¥h WjWk, 
d, Xj/j and xjfj, and remaining non-derivative terms. By equating the coefficients of the 
above terms to zero, one obtains the following system of coupled differential equations 

dC ^ ( dC' dD' \ ^ 

C F + 2a^F + a 2 F ; T -^ £i¥l + —i £i \j/* + aF ; y(C i £ i x l / I +D I £ i x l /*)=0, 
da \da da J r-[ 

(32) 





r\w^> + w™)=°' F iw^ + wr^ - <35) 



4 dd 
3Coy/j + aDj + a ^ I 

' 1 \ J 



dDi 



0, 



(37) 



3CbyJ +aCj-a^ 



dDi „ 



0, 



(38) 



4 

3C V + aV' £ (Qe; ^ + y/f ) = 0. (39) 

!=1 

In equations (1321) through (139b it was introduced the symbol 

£,- = +1 for i = 1,2; 

£j = -l for i = 3,4. (40) 

The system (I32l)-(l39l) is solved for two cases: F' = and F' ^ - the details of 
the integration of this system can be found in the reference [|3~2ll . The corresponding 
solutions are exposed as follows. 

Case F' = 0. 

F = constant, V = X X ¥, (41) 
with k, /3 , A= constant. 
Case F' ^ 0. 



Set 1 



3 3 
C = ka, Cj = ~2 k ¥j + P £ j¥j, D J = ~2 kVj ~ $ E Wh 

F = a x ¥, V = AY, (42) 

Set 2 



C = -, Cj = -^Vj+Pejifi, Dj = --ty-peM, 

F = a x ¥ l/ \ V = A V I', (43) 



where a is a constant. 



Solving the field equations 

We cannot distinguish physically the potential V given by ( l4Tjb from a mass term 
- see (PT3l . Then one can consider V = AY = mW, so that A corresponds to a mass 
parameter. Thus from now on A will be replaced by m. 



From (1241) and (1231 ) one can build an evolution equation for X V, which reads 
¥ + 3/^ = 0, so that 



(44) 



where is a constant. Then, once *F = ^(a), equation (1291 becomes a function only of 
a and, as it will be seen below, it can be directly integrated for all the cases. Therefore it 
is not necessary to calculate the constants of motion. 

Case F' = 0. The energy density and pressure of the field follow from (l27l and (1281) 
through l44l yielding 

P/=-^r> Pf = °> ^ 

which imposes that > 0. 

The choice F = constant =1/2 corresponds to a minimal coupling between the spinor 
and gravitational fields. In so doing, by inserting the obtained pf into (12~9~1 ), it furnishes 
the scale factor 



2 ) ("J/ 



(46) 



which describes a decelerated Universe as it was dominated by a standard pressureless 
matter field. 

Substituting (1461) in (1241) . we have the solution for the spinor field 



V(0 = a 



3/2 



\j/^ e -" nt 
yO e imt 

\ y%e imt ) 



,-t 



(47) 



with the y/J being constants such that = Vi° ¥i + V4°V^ ~~ Vf ^ ~~ wfw®- 

Case F' 7^ 0. For the Set 1 solutions, when F = a^, the energy density and pressure 
read 

P/ = "2^' ^/ = "P/' (48) 



so the integration of (1291) gives 



a(?) 



„Kt 



where 



K 



m 



12a' 



(49) 



which trough (1241 renders the following solution 

Xi„—2imt 



v(0 



yQ e -2imt 



2unt 
,2imt 



3K. 

e 



(50) 



with the y/j following the same constraint as for solution (l47l ). 

We infer from (l48T h that one must have < 0, which imposes that a < since the 
coupling F = a v I / is a positive quantity. And by its turn, the condition a < implies that 
K ia a real quantity. So the solution (|49l ) can describe a De Sitter-like Universe. Hence 
in this case we could identify the spinor field with the "inflaton". The result (148T b asserts 
that the pressure is always the negative of its corresponding energy density, similarly to 
the cosmological constant state equation. 

For the Set 2 solutions, when F = am 1 / 3 , equation (|29l) does not have a solution. 



SCALAR AND SPINOR FIELDS MINIMALLY COUPLED WITH 

GRAVITY 



Consider a general action for a scalar field and a spinor field minimally coupled with 
the gravity 



where stands for the Lagrangian of a common matter field. 

For a flat F-R-W metric applied to (ISTI) . with the fields and i// being spatially homo- 
geneous and the common matter pressureless, we can write the point-like Lagrangian 



where denotes the energy density of the common matter at an initial instant. 

Imposing that the energy function associated with (l52l) is null, the result is the Fried- 
mann equation 



General action and field equations 




(51) 



Jz? = 3ad 2 -a 3 



Uj + l -a 3 (Wfy-xj/fy) +a 3 V + p° M , (52) 




(53) 



with py and denoting the energy densities, which are defined as 



1 . 7 

Py = V, p0 = -r + t/. 



(54) 



From the E-L equations for a applied to (1521) . we have the acceleration equation 

a_ pM + P V + P<l>+3(P\ir + P(t>) 



(55) 



and for i// and i//, the Dirac equations for the spinor field and its adjoint coupled with the 
gravitational field, respectively, 

3 n dV 

Hr+-H\if + if\if— = 1 (56) 

y+-Hy-iyf— = 0. (57) 

In (1551) the pressures are defined as 

dV S 2 
PY = *jy- V > Pt = Y~ U - (58) 

The E-L equations for furnishes 

. dU 

$ + 3H$ + — = 0, (59) 
d(p 

and so we have the Klein-Gordon equation for the field 0. 



Determination of the Noether potentials 

Expressing the Lagrangian (1521 in terms of the spinor components, it takes the form 

if = 3ad 2 - a 3 (^--U^j+ l -c? £ i# - yif ft) + a 3 V + Pm- (60) 

Thus the configuration space of the system is represented by (a, 0, y/j, i/z^d, 0, y/^, y/y). 

The Noether symmetry condition, L^S£ = XJz? = 0, applied to (l60l) . with X defined 
for our problem as 

d . d d . d * / d d . d . d \ 

da da d(j) d<j> ; t"J \ dx/f* d\j/j d\j/* a%) 

renders the following system of coupled differential equations 

C + 2a^r— = 0, 3C + 2a^r— =0, 6^— -a ^— = 0, (62) 

* / dd dDi \ ^ /dC ; - dDi A n 



z=l 



<9a da ' J ^1 V ^0 ^0 



A f dCi dDi \ n 
3C Q¥j + aDj + a^ = 0, 



(64) 



3C \if]+aC j -a^l^W--^Vt) =0, (65) 

= ° 5 3 — = 0, 3— ^ = 0, - — = 0, (66) 

^v; 3vo ^v? ^ 

3C (£7 + V) + aD ^- + a Y (Qe^ + D^-yf ) ^ = 0. (67) 



One can see from equations (1661 ) that the coefficients Co and Do are functions only of 
a and 0. Then, assuming that Co and Do are separable functions 

C = ci(a)c 2 ((j)), D = d l (a)d 2 ((j)), (68) 
we obtain the solution for the system (I62l)-(l67l) 



Ae a *+Be- a * 

C = f= , (69) 

V6a 



Ae a *-Be- a + 
D o = ^372 > ™ 

Cj = -c ^ Perf, (71) 
Ae^ + Be~ ai> 

D J = ~ a ^372 V0' - PWh ( 72 ) 

U = Uo(Ae a $ -Be- a ^ 2 , V = V x lf, (73) 

with A, B, Uo, Vq and /3 being constants and a = \ \J\. 

Dark matter as femionic particles. The Noether potential V = Vo^ is essentially a 
term of mass. Then, as before, we replace the constant Vq by m. And from (j54j) i and 
(j58|) 1 the energy density and pressure are 

p v = rrW, p v = 0, (74) 

characterizing a field of pressureless matter - when *P is a non-negative function. 

This result suggests that the spinor field behaves as a standard matter field. But 
obviously this field has a nature different from that of the common matter since it 
describes a "fluid" exclusively composed by fermionic particles. Such a field produces 
an additional pressureless matter that may be identified with the dark matter. For this 
propose, we must assume that these fermionic particles interact only gravitationally 



with the common matter or have very weak non- gravitational interactions with it. On 
the other hand, the potential (173]) i can produce an accelerated expansion. Therefore we 
may identify the field with the dark energy. Hence, from this point on, the dark sector 
will be identified with the fields and y. 



Analytical solutions of the field equations 



Taking U and V given by (1731) for (I53l)-(159T), we have the following system to solve 
3H 2 = ^-+m^ + ^- + U (Ae a $ - Be~ a *} 2 , (75) 

H 2 + 2- + - Uo( Ae a * -Be- a *) 2 = 0, (76) 

a 2 V / 

$ + 3H<j) + 2aU (A 2 e 2a ^ - B 2 e- 2a ^ = 0, (77) 

3 3 
\j/+-H\j/ + im'f\j/ = 0, 1j/ + -Hy - imxj/'f = 0. (78) 

Once the dynamical system now presents a Noether symmetry, one has an additional 
equation provided by the constant of motion, namely 

M » = C "^ +D "W + ^X Cl Bwf +Di S*) 

= (Ae a * +Be- a ^ Vfad+ (Ae"* - &T°*) a 3 / 2 <p + a 3v P, (79) 

which is determined from ©. In this equation one can make /3 = ifio, with j8q being a 
real constant, in order to have a constant of motion whose value is real. 

Consider now a transformation of variables that changes the set of coordinates 
{a, 0, y/j, \ffj} to {z,u,vj,Wj}, which satisfies the system 

du du Jt, ( du du \ 

c »5S +D »50 + S( Ci aW +D '^) =0 ' <80) 

with z being the cyclic variable. This system is obtained by applying ©. 



Using the constant of motion (|79l ) and a transformation of variables satisfying (f80l) - 
(1831 ), we will look for a solution to the field equations (|75l)-(|78T). 

Firstly, equations (|78|) i and (TTSlb can be reduced to an equation for *F with the same 
form of (|44)) . then the form of the bilinear is again ¥ = *I> / a 3 , so that ¥ > by (HU). 
And by solving (|78]) i in terms of a, the solution for \y reads 



a~ 3/2 , (84) 



2 



with the subjected to the same constraint of (I47T ). The solution *F = ^o/^ 3 through 

(|74|) implies the energy density = mM'o/a 3 , so the extra matter field is naturally 
added with the common matter field, which corroborates the possible identification of 
the spinor field with the dark matter. 

Since equations (1781) can be independently solved in terms of a, we effectively have 
a reduction of the system (T75"1)-(|78T) to equations that involve only the variables a and 0. 
Thus we will solve the problem through the reduced system 

3H 2 = P°M+fQ + f_ + [Jo U e a, _ Be -a,\\ (g5) 
2 V / 

(p+3H<j>+2aU (A 2 e 2a ^ -B 2 e- 2a ^ = 0, (86) 

(Ae"* +Be- a ^ V&ia+ Ue a * -Be~ a ^ a 3 l 2 § = M , (87) 

where M =M + jSo^o and the result = x ¥ /a 3 was used. 

Hence we need to find a new set of variables related only to a and 0. An adequate 
particular solution of this type for the system (l8~0l-(l8~3T) is 



u = a 3/2 = , v, = w: = 0, z = a 3/2 = . (88) 

V6AB 3 J V6AB 

Using (1881) we can express the constant of motion in terms of the cyclic variable 
. M 



4AB 1 



so that z(t) = zit + Z2, (89) 



with z\ = Mq/4AB and zi = constant. 

Equation (1751) takes the following form in the new variables 



z 2 = u 2 + 3ABU u 2 + ^±^. (90) 

TAB 

Substituting z from (l8~9l) in the above equation, we get its integration and obtain 



u(t) = uosm((Ot + bo), 



(91) 



where 



M 2 - SAB (pQ +m^ ) 
48A 3 5 3 [/ 



Uq 



(0 = 3ABU , 



(92) 



with bo being a constant. 

From the solutions z(t) and w(?), expressed in terms of the original variables through 
the relations (1881) . we get the explicit forms of a(t) and <j>(t) 

a(t)= (J^J ' {z\t 2 + 2ziz 2 t+zl-ulsm 2 {m + b Q )} ll \ (93) 



1 ( z\t + Z2 + u sin (cot + b ) \ 1 



In 



hi 



2a [zit + Z2 — Mosin(ct)f + bo) J 2a 
Now from (|84l ) one can write the time evolution of the spinor field 



(94) 



¥(t) 



V2U 

CO 



y®e" nt 
\ \$e mt 



\ 



J 



(95) 



-1/2 



where 

©(t) = {z 2 t 2 + 2ziz 2 t + zl-ulsm 2 (cot + bo)y (96) 

The solution (|93l describes a Universe with an oscillatory expansion rate and can 
account for the transition from a decelerated to an accelerated era, in according to the 
observations. Further, it delineates a Universe that will return to a decelerated period 
and turn to accelerate in the future and so on, such that these alternated periods of 
acceleration and deceleration have a behavior like that from a deadened oscillator - when 
t — > oo, one has a — > 0. The time evolution of the scalar field (l94l) makes the pressure 



(|58|)7 oscillate between a negative and a positive value, and so the field alternately 
behaves as dark energy and matter, being the cause of the oscillatory expansion. The 
detailed analysis of the above solutions can be found in the reference Q33D - 



CONCLUSIONS 

Generic cosmological models with spinor and scalar fields were constrained through 
the Noether symmetry method. The proprieties of the Noether symmetry rendered the 
complete integration of the dynamical systems. The cosmological solutions showed that 
the spinor field can play the role of "inflaton" - non-minimally coupled case -, such that 
its equation of state is similar to that of the cosmological constant; and for the minimally 
coupled case, a behavior of standard matter was produced. For the model with spinor and 
scalar fields minimally coupled, the spinor field presented a standard matter behavior - 
which may describe the dark matter - and the scalar field an oscillating equation of state, 
which passed from a matter-like to a dark energy-like description producing alternated 
periods of accelerated and decelerated expansion. 
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